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In this short review we summarize experimental information and theoretical results for the low- 
energy dissociation cross sections of charmonia by light hadrons. These cross sections are required 
for the simulation of charmonium absorption through collisions with comovers in heavy ion collisions, 
which competes with quark-gluon plasma production as a charmonium-suppression mechanism. If 
the cross sections are sufficiently large these dissociation reactions may be misinterpreted as an 
effect of quark-gluon plasma production. Theoretical predictions for these RHIC-related processes 
have used various methods, including a color-dipole scattering model, meson exchange models, 
constituent interchange models and QCD sum rules. As the results have been largely unconstrained 
by experiment, some of the predictions differ by orders of magnitude, notably in the near-threshold 
regime that is most relevant to QGP searches. 

PACS numbers: 12.38.Mh, 13.75.-n, 14.40.Lb, 25.75.-q 



INTRODUCTION 

Many unusual subjects have been studied in the name 
of QCD. One of the more unusual, which has arisen in 
the field of heavy ion collisions, is the size of cross sec- 
tions of charmonia on light hadrons. This has attracted 
attention because of the proposal by Matsui and Satz Q] 
that suppression of J '/tp production could be used as a 
signature for the formation of a quark-gluon plasma. 

This suggestion, like many signatures proposed for the 
quark-gluon plasma, is perhaps excessively intuitive. The 
idea is that a QGP will screen the linear confining inter- 




FIG. 1: We wish to distinguish between 1) weak and 2) strong 
J/ip absorption by light comovers. 



action between quarks, so that a cc pair produced within 
a QGP will be less likely to form a bound cc charmo- 
nium resonance, as in Fig.l, but instead will more likely 
separate to form open-charm mesons. 

Even if this simple picture of cc production in a QGP 
is qualitatively correct, it can only be confirmed easily if 
the competing direct charm production and scattering by 
the initial relativistic nucleons is understood and 
if there is little subsequent dissociation of the charmonia 
by the many other "comoving" light hadrons produced 
in such a collision. To summarize the last point, if char- 
monium + light hadron "comover" dissociation cross sec- 
tions are small (case 1, top of Fig.l) and the background 
of direct charm production from the initial nucleons is 
understood, one may have a useful signature of QGP 
formation, but if the comover dissociation cross sections 
are large (case 2, bottom of Fig.l) one must distinguish 
a QGP-reduced charmonium production amplitude from 
subsequent dissociative scattering, and the interpretation 
of this signal will therefore be ambiguous. 

Thus it is of great relevance to the interpretation of 
RHIC physics to establish the approximate size of these 
low-energy cc + light hadron cross sections. 

EXPERIMENT, OR WHAT PASSES FOR IT 

Unfortunately we have no charmonium beams or tar- 
gets, so the experimental cross sections must be inferred 
indirectly and are poorly known. The earliest estimates 
of lower energy charmonium hadronic cross sections came 
from J/vb photoproduction experiments in the mid 1970s, 
which were interpreted in terms of a J/ip+N cross section 
given additional theoretical assumptions. Early Fermilab 
and SLAC photoproduction experiments gave rough esti- 
mates of ~ 1 mb for ctjm+Nj assuming vector dominance, 
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FIG. 2: A fit of Eq.(l) to experimental J/ip production versus 
path length [IJ; the line corresponds to 6.2 mb. 



for photon energies from E 1 « 13 to 200 GeV HQ. A 
subsequent SLAC photoproduction experiment in 1977 
used the A dependence of J/ip absorption to estimate 
a rather larger cross section of ctj/^+n = 3.5(0.8) mb 



at £ 7 k 17 GeV (y/s w 6 GeV) [jj. The vector domi- 
nance hypothesis may have lead to an underestimate of 
the cross section in the earlier references 

In heavy ion collisions these cross sections may be es- 
timated from the ratio of lepton pairs produced in the 
J I if) peak to "background" Drell-Yan pairs nearby in en- 
ergy. Since the J /if) must reach the exterior of the nuclear 
target to decay into a sharp mass peak, this ratio gives 
us an estimate of the absorption cross section through a 
classical survival probability formula, 

<j{J/if> — > A* M~ )/c(Drell — Yan p, + p~) 



exp(-p a a j% +N L) 



(1) 



where p is the mean nucleon density and L is the es- 
timated mean path length in the experimental nuclear 
system. A "naive" interpretation of the J/if> production 
data from collisions of various nuclear species using this 
formula gives er°fe +N ~ 6 mb at *Js ~ 10 GeV Q, with 
a numerically similar result for the if)' . 

Of course one may raise many questions about the va- 
lidity of this simple estimate, including the importance 
of shadowing in Drell-Yan, the neglect of J /if) scattering 
by other light hadrons formed in the collision (for exam- 
ple 7r and p), and the assumption of a single, constant 
J/i/j+N cross section in all circumstances. 

In addition to these more or less direct measurements 
there have been several attempts to infer charmonium 
cross sections from related processes, given additional 
theoretical assumptions. A series of theoretical papers 



[ToL Ull ll2L Il3| has estimated the rather weak closed- 
flavor cross sections such as J/ip+w — > ifi'+ir near thresh- 
old from the observed dipion decays (here if)' — ► mrJ/ip). 
Redlich et al [l4| have used charm photoproduction data 
combined with a vector dominance model to estimate 
the J/ifM cross section, and find quite small values near 
threshold. Hufner et al. 0] again argue however that 
vector dominance is not justified for this process, so that 
the J /if) + N cross sections estimated in this manner are 
inaccurate. 

Recently, concerns have been expressed that the J /if) 
and if)' wavefunctions have not had sufficient time to 
form within the nucleus in these collisions, so experi- 
ment may instead be measuring the cross section for a 
small initial cc "premeson" on a nucleon. One can in- 
crease the time spent in the interior of the nuclear sys- 
tem by selecting small and even negative xp events, as 
has been done by E866 at Fermilab. As discussed by He, 
Hufner and Kopeliovich [lil Il7| , this leads one to infer 
cr}^ +N = 2.8(0.3) mb and a$' +i , = 10.5(3.6) mb respec- 
tively, also at y/s » 10 GeV. This is rather more satisfy- 
ing to people who have an intuitive notion that the larger 
if)' should have a larger reaction cross section. Actually 
the connection between cross section and physical extent 
is less direct (compare KN and KN), and in any case 
the relative proximity of inelastic thresholds alone would 
suggest a larger if)' cross section. These experiments also 
indicate a preference for dissociation over elastic cross 
sections in this energy region by roughly a factor of 30 
00. 



THEORY: INTRODUCTION 

To quote B.Muller in 1999, "...the state of the theory 
of interactions between J/ip and light hadrons is embar- 
rassing [...]. Only three serious calculations exist (after 
more than 10 years of intense discussions about this is- 
sue!) and their results differ by at least two orders of 
magnitude in the relevant energy range [...]. There is 
a lot to do for those who would like to make a serious 
contribution to an important topic." |l8j . 

The theoretical situation has improved considerably in 
recent years, at least in terms of the number of calcula- 
tions if not in the understanding of the scattering mecha- 
nism. A list of cc + light hadron cross section calculations 
is given in Table I. 



Color Dipole Model 

The most cited work, albeit furthest in its predictions 
from a low-energy "theoretical mean" , is the color-dipole 
scattering calculation of Kharzeev and Satz 19]. This 
reference is basically a restatement of the color-dipole 
model developed in the late 1970s by Peskin and Bhanot 
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Method 


Init. State 


Ref. 


color dipole 


j/V>+n 


02 




J/i> + (tt,N) 


[20] 






[15, 21J 


meson ex. 


J /if) + IX 


[27] 




J/lp + TT yP 


[28, 22, 30, 31, 32] 




J/V'+N 


[33, 34] 




J/lp + 7T, K, p,N 


[35, 36J 


constit. int. 




[27, 37] 






[22] 






[22J 




(J/ip,i>') + (n,p) ; XJ + tt 


[40] 




J/V-+N 


[4JJ 




(JM^) + (7T,N ) 


[42] 


QCD sum rules 


j/y> + 7T 


[43, 44] 



TABLE I: A summary of near-threshold cc + light hadron 
cross section calculations. 



for scattering of light hadrons by Coulombic 
bound states of very massive quarks. According to Pe- 
skin, the criterion for validity of this approach is "...not 
met even for the bb system." |25|. so there may be large 
systematic errors at the cc mass scale. This approach cer- 
tainly makes marginal approximations for charmonium, 
such as the use of Coulombic wavefunctions (which are 
far from accurate for cc) and the introduction of a QQ 
binding energy (which is hard to interpret for charmo- 
nium, and is taken to be 2 Md — Mj/^ by Kharzeev and 
Satz). These color-dipole scattering formulas also im- 
plicitly assume that charmonia are small relative to the 
natural QCD length scale. Since potential models actu- 
ally find rms cc separations of about 0.4 fm for the J/ip, 
0.6 fm for the Xc states and 0.8 fm for the ip' [26j, this is 
also a dubious assumption. 

Although this approach has problems with justifica- 
tion for cc, the predictions are nonetheless interesting as 
estimates of the scale of these cross sections assuming 
a color-dipole scattering mechanism, and the approach 
could presumably be extended to lower energies by gen- 
eralizing the wavefunctions and interaction. The formula 
for the J/ifj+N cross section quoted by Kharzeev and 
Satz is 

<jj^ +N = 2.5 mb - ^\ (2) 

where A = (s — Mj,^ — M^)/2Mj/^, the constant A 
is defined to be " ~ Mn + eo" according to the text 
following Eq.(24) of |HJ (we assume the equality), and 
the "binding energy" eo is set to 2 M-q — Mji^,. The result 
is shown in Fig. 3, together with the single lower-energy 
SLAC experimental point Q. 

Evidently the Kharzeev-Satz cross section is smaller 
than this SLAC point (which was an inferred cross sec- 
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FIG. 3: The Kharzeev-Satz J/ip+~N total cross section and 
the 1977 SLAC result 0. 

tion and certainly needs confirmation) by about two or- 
ders of magnitude, and falls precipitously as ^/s is de- 
creased. Their calculation actually has no direct in- 
formation about physical thresholds, so Kharzeev et al. 
typically leave their curves "dangling" just below s/s — 
5 GeV. (See Fig. 2 of 0] for example.) If we plot their 
formula Eq.(2) for J/^+N at low energy (Fig. 3), we find 
the unphysical prediction of a zero cross section near 
4.5 GeV, whereas the physical A c + D threshold is at 
4.15 GeV. Obviously this calculation is inapplicable at 
low energies, which is unfortunate because this is the 
regime of greatest interest for QGP searches. 

Similar cross section calculations have since been re- 
ported by Arleo et al. [2(J and Oh et al. _2lj (also 
using Coulomb wavefunctions, but incorporating physi- 
cal kinematical constraints) and by Hiifner et al. [l5j, 
using a similar dipole scattering model with more re- 
alistic cc wavefunctions. Oh et al. |2l| note that this 
model predicts an extremely large cross section ratio of 
a^+p/aj/^+p - 2000-5000 at 4.2 GeV. Hiifner et al. find 
J/ijj + p and ip' + p cross sections of 3 mb and 10 mb 
respectively at -^/s = 10 GeV (see their Fig. 10), and do 
not consider use of this approach justified near threshold. 

Meson Exchange Models 

Several calculations of charmonium + light hadron 
cross sections have been reported assuming i-channel 
charmed meson exchange. Of course this picture is also 
problematic, since the range of the exchanged charmed 
meson would be only about 1/Md ~ 0.1 fm, and the as- 
sumption of nonoverlapping hadrons at this separation 
is clearly invalid. (This is the Isgur-Maltman ^5| ar- 
gument as to why vector meson exchange is unjustified 
as the source of the short-ranged NN core interaction.) 
Nonetheless it is again interesting to see what scale of 
cross section is predicted by this type of model, since it 
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FIG. 4: The Matinian-Miiller t-channel meson exchange re- 
sults for J/ip + tv and J/tp + p inelastic cross sections 28]. 



might at least incorporate the correct scales and degrees 
of freedom, and it assumes a different scattering mecha- 
nism from the color-dipole model advocated by Kharzeev 
and Satz. 

The first such meson exchange calculation, due to Ma- 
tinian and Miiller [2^, considered i-channel D exchange 
as the mechanism for the reactions J/^p + n — > D*D + h.c. 
and J ftp + p — > DD; their results are shown in Fig. 4. 
Note that 500 MeV above threshold these cross sections 
lie in the 0.5 to 1 mb range. A subsequent calculation by 
Haglin _35], who assumed an SU(4) invariant gauge field 
effective meson lagrangian, found few-mb scale results for 
these cross sections near threshold. 

Similar meson exchange dissociation cross sections cal- 
culations for J/i/j+N were reported by Haglin j^] , Sibirt- 
sev, Tsushima and Thomas and Liu, Ko and Lin [34| . 
Haglin found a peak of about 7 mb near y/s = 4.3 GeV, 
whereas Sibirtsev et al. found a peak cross section of 
about 1 mb near y/s = 4.6 GeV (Fig. 5). Liu et al. found 
a scale similar to Sibirtsev et al. near threshold, but 
concluded that production of charmed meson pairs dom- 
inated at higher energies. Since Kharzeev and Satz found 
that the color-dipole scattering model predicts negligible 
J/tp+N cross sections at low energies (Fig. 3), evidently 
flavor-exchange processes are predicted to be dominant 
near threshold. 

Subsequent meson-exchange studies by Lin and Ko 
[H, Haglin and Gale [l^ and Oh et al. [Mill found 
that the assumption of pointlike mesons in the effec- 
tive lagrangian leads to unrealistically large cross sections 
even at moderate energies, so that hadronic form factors 
must be incorporated in the calculations. With plausi- 
ble but rather arbitrary form factors these cross sections 
are greatly reduced, so that they are again found to be 
typically of few-mb scale (see Fig. 6 for example). 




(GeV) 



FIG. 5: The t-channel meson exchange cross sections for 
J I ib + N — > A c + D found by Sibirtsev, Tsushima and Thomas 
|33|. The smaller contribution is from D exchange and the 
larger is from (non-interfering) D* exchange. 
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FIG. 6: A strong suppression of dissociation cross sections 
is found on incorporating hadronic form factors in meson ex- 
change models. This example is Fig.4 of Lin and Ko |2fll |. 



In summary, the current, rather unsatisfying state of 
affairs in meson exchange models of charmonium dis- 
sociation is that mb-scale cross sections are anticipated 
near threshold, but their precise values depend on poorly 
known hadronic form factors, as well as on effective me- 
son lagrangians of unknown accuracy. Future studies 
that can put the effective lagrangians and form factors 
on a more sound theoretical basis would be an important 
contribution to this approach. One interesting possibil- 
ity, which has recently been investigated by Navarra et al. 
and Matheus et al., is to derive the hadronic form factors 
from QCD sum rules [iil . An especially attractive 
possibility explored very recently by Deandrea, Nardulli 
and Polosa |4jjj is that the hadronic form factors may be 
evaluated explicitly in terms of quark model wavefunc- 
tions. 
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Constituent Interchange Model 

In this approach one uses explicit nonrelativistic quark 
model wavefunctions for the external hadrons and calcu- 
lates the cross section assuming a constituent interchange 
scattering mechanism, driven by the Born-order matrix 
element of the standard quark model Hamiltonian. Con- 
stituent interchange cc + qq cross sections have their 
strongest support just a few hundred MeV in y/s above 
threshold, since the overlap integrals are damped by the 
external meson wavefunctions at higher momenta. 

This technique, which has no free parameters once 
quark model wavefunctions and the interquark Hamil- 
tonian are specified, has been shown to compare reason- 
ably well with experimental low-energy hadron-hadron 
scattering data near threshold for a wide range of 
annihilation-free reactions. There are four quark inter- 
change diagrams in meson- meson scattering (Fig. 7), each 
of which has an associated overlap integral of the exter- 
nal meson wavefunctions convolved with the interquark 
Hamiltonian. Constituent interchange is forced at Born- 
order because Hi a A° ■ A a changes each initial color- 
singlet qq meson into a color octet, which has overlap 
with color-singlet final state mesons only after quark 
line interchange. Feynman rules for these diagrams were 
given by Barnes and Swanson |4flj . 
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FIG. 9: tp' + ty cross sections |4(j. Note the change of scale 
relative to Fig. 8. 



T1 T2 

FIG. 7: The four constituent interchange scattering diagrams 
evaluated in cc + qq cross section calculations (prior for- 
malism). The "exchange" is the full quark-quark interaction 
Hamiltonian Hi. 

The first charmonium cross section calculation using 
this approach was due to Martins, Blaschke and Quack 
[37I. who considered the reactions J/ip + tt — > D*D + h.c. 
and D*D* (The amplitude for J/t/j + tt — > DD is zero in 
the nonrelativistic quark model without spin-orbit forces, 
and has been found to be quite weak in a relativized cal- 
culation (2t| •) Martins et al. found that these exclusive 
final states have numerically rather similar cross sections 
(except for their different thresholds), and give a maxi- 
mum total cross section of about 7 mb at yfs 4.1 GeV. 
Wong et al. have carried out similar constituent inter- 
change cross section calculations [3^, l^, ^3 1 using nu- 
merically determined Coulomb plus linear plus smeared 



hyperfine quark potential model wavefunctions, with pa- 
rameters fitted to the qq and cc meson spectra. Figs. 8-10 
show some recent results from Ref. • The J/ip+ir cross 
section is somewhat smaller than was found by Blaschke 
et al., and peaks at about 1.4 mb just above 3.9 GeV. 

The difference lies mainly in the treatment of the con- 
fining interaction; Blaschke et al. treated confinement as 
a color-independent Gaussian potential between quark 
and antiquark (hence they include only diagrams CI and 
C2), whereas Wong et al. used the conventional A a • A a 
form between all pairs of constituents. This leads to de- 
structive interference between the C and T diagrams due 
to opposite color factors, and hence to a much reduced 
total cross section. Ref. |4C| also considers tp' + n scatter- 
ing; the rather large cross section found for this process 
is shown in Fig. 9. As a final example, the constituent- 
interchange J/tp + p cross section from Ref. 0] is shown 
in Fig. 10; this diverges at threshold for the simple reason 
that the DD channel is exothermic. 
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FIG. 10: J/ip + p cross sections |4C 
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FIG. 11: QCD sum rule results found by Navarra et al. [IIJ 
for the cross sections J/ip + it — > DD, D*D + h.c. and D*D* 
(distinguished by the thresholds) and their sum (solid) . 



It is interesting that the simple two-parameter function 

<r(s) = o max {e/e max f e KW^ ) , (3) 

often provides an accurate parametrization of the con- 
stituent interchange cross sections. In this formula e = 
y/s — Mc — Md is the cm. energy above threshold and 
Vmax is the cross section maximum, at e max . The thresh- 
old exponent p is determined by the angular quantum 
numbers of the hadrons, and is ±1/2 + L^ n (for en- 
dothermic/exothermic), where L^ n is the lowest angular 
momentum allowed for the final meson pair CD. 

The much more complicated problem of charmonium- 
nucleon scattering has also been investigated in the con- 
stituent interchange model. J/tp+N and ^i'+N cross 
sections were derived by Martins |42j using a simplified 
quark-fdiquark model of the nucleon; this approach gave 
several-mb peak cross sections not far above threshold. 
Black |4l| nas carried out the full J/ip+p constituent 
interchange calculation for several final states, assum- 
ing hyperfine, Coulomb and linear interactions. He finds 
that the final state D°A+ is dominant, with a surprisingly 
large (ca. 12 mb) peak cross section just above threshold. 

QCD Sum Rules 

The application of QCD sum rules to the determina- 
tion of charmonium dissociation cross sections, due to 
Navarra et al. ^| and to Duraes et al. [SJj is a rela- 
tively recent development. This method relates the scat- 
tering amplitude to a sum of operator vacuum expected 
values, and gives a model independent result to the ex- 
tent that these expected values are known experimentally 
and the set of operators chosen does indeed dominate the 
scattering amplitude over the chosen kinematic regime. 
There are also systematic uncertainties in the approach 
due to the details of a Borel summation and treatment 
of a continuum contribution to the amplitudes. 



The studies published to date are a calculation of the 
J/ij) dissociation reactions J/ip + n — > DD,D*D + h.c. 
and D*D* 0], followed by a more detailed investigation 
of the same processes The sum rule results (Fig. 11) 
appear to confirm the mb-scale of near-threshold dissoci- 
ation cross sections, in qualitative agreement with both 
meson exchange and constituent interchange models. 



CONCLUSIONS 

We have reviewed theoretical predictions and related 
experimental results for the low-energy dissociation cross 
sections of charmonia on light hadrons, which are of 
great importance for QGP searches in heavy ion colli- 
sions. Four theoretical approaches have recently been 
applied to this problem, which are a color-dipole scat- 
tering model, meson exchange, constituent interchange, 
and QCD sum rules. Near threshold the color-dipole 
model predicts very small cross sections, whereas meson 
exchange, constituent interchange and QCD sum rules 
all predict mb-scale cross sections. At present there are 
no direct measurements of these dissociation cross sec- 
tions near threshold. It would clearly be of great inter- 
est to measure any of these low-energy dissociation cross 
sections experimentally, both for the relevance to QGP 
searches and as a valuable test of the theoretical mod- 
els of hadron scattering that have been applied to this 
problem. 
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